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This letter calculates the contribution of electron-phonon interaction to thermoelectric effects in 
graphene. One considers the case of free standing graphene taking into account interaction with 
intrinsic acoustic phonons. The temperatures considered range from liquid nitrogen to the room 
level. It has turned out that the contribution to thermoelectromotive force due to electron drag by 
phonons is determined by the Fermi energy in the sample and phonon relaxation time. The explicit 
temperature dependence of the contribution to thermoelectromotive force deriving from by phonons 
is weak in contrast to that due to diffusion, which is directly proportional to temperature. The 
dependence obtained suggests that at the temperatures considered, a high carrier concentration and 
for samples with a high thermal conductivity, the phonon drag contribution can become dominant. 
Thus a theoretical limit has been established for a possible increase of the thermoelectromotive force 
through electron drag by the intrinsic phonons of graphene. 

PACS numbers: 72.80.Vp, 63.22.Rc, 81.05.ue, 73.50.Lw 


The transport properties of graphene have been at¬ 
tracting considerable interest in recent years^P— . A 
study was performed of the electrical^i^ and thermal 
conductivity2r2£, thermomagnetic effectsii and electrical 
phenomena generated by irradiatio n^^d^ . 

Particularly promising for the field of applications ap¬ 
pear to be experimentali^^— and theoreticali2r.2i stud¬ 
ies of the thermoelectric properties of graphene. The 
high electrical conductivity combined with a significant 
Seebeck coefficient (thermoelectromotive force or ther¬ 
mopower), provide a basis for expecting graphene to re¬ 
veal a fairly high thermoelectric figure of merit. Scatter¬ 
ing from phonons is essential for evolvement of electron 
transport in graphene^i^i^^— . Thus, in graphene, as in 
other sp^ carbon materials^^— , one can envisage ampli¬ 
fication of thermoelectric effects as a result of the phonon 
drag effect discovered in 1946 by L.E. Gurevich^i. 

We calculate the electric current generated by phonon 
drag in the presence of a temperature gradient in a 
graphene sample, obtain the corresponding contribution 
to the Seebeck coefficient and compare it with the con¬ 
tribution of diffusion to thermopower. 

The conditions chosen for the calculations are the most 
favorable for experimental observation of the effect of 
electron drag by phonons. The assumptions accepted 
for the sample and external conditions in the system are 
as follows. We consider electrons near the Dirac point 
having an idealized linear dispersion law = Hvpk (the 
Eermi velocity vp ~ 10®cms“^). The Eermi level ep in 
the system is assumed to be much higher than the tem¬ 
perature, opening a way to application of the degener¬ 
ate electron gas approximation, which implies monopo¬ 
lar conduction. The condition of the onset of degener¬ 
acy defines the upper limit of temperature below which 
this theory is applicable. The absolute value of the elec¬ 
tron wave vector k near the Fermi level will be denoted 


in what follows by kp. We are going to consider a 
free-standing graphene monolayer and take into account 
the contribution of the intrinsic acoustic in-plane (LA 
and TA) phonons only, which correlates with previous 
studies^ of the mechanisms underlying electron relax¬ 
ation in graphene. The dispersion law of such lattice 
vibrations is given by fiWph{q) = hqvg, where sound ve¬ 
locity is ^ 2 • 10®cm s“^. 

On the other hand, it is assumed that the temper¬ 
ature is high enough to maintain an equipartition (EP) 
regime of electron-phonon interaction, in which the num¬ 
bers of occupation of the phonon modes contributing to 
the drag effect are larger than 1. The phonons interact¬ 
ing with electrons have a wave vector on the order of kp 
and energy hvskp. 

In the case of high temperatures, the equilibrium dis¬ 
tribution function for phonon mode occupation numbers 
transforms into (q) = T/Hu!ph{q) for both phonon 
emission and absorption processes. This defines the lower 
temperature limit of applicability of the theory devel¬ 
oped. Only the contribution from acoustic phonons is 
considered, and, hence, the occupation numbers of op¬ 
tical phonons must be small. This condition specifies 
the second upper limit in temperature for the theory 
developed. Fig. 1 demonstrates all characteristic ener¬ 
gies, wave vectors and velocities related to the problem 
being analyzed. We readily see that in actual fact the 
above conditions are satisfied by the very broad temper¬ 
ature range, including room temperature. Finally, we 
neglect the reverse influence of the electron subsystem of 
the crystal on its phonon subsystem. 

The Fermi energy of graphene ep is known to be 
related^ with carrier concentration rig as ep = hvpffpng. 
In transistor-based electronic devices the concentration is 
controlled by properly varying the gate voltage and is de¬ 
termined by the capacity of the structure. These param- 
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FIG. 1. Schematic of the energy scales essential for the prob¬ 
lem under study. Dispersion of electrons (black line) and dis¬ 
persion of acoustic phonons (grey line) are shown on an ar¬ 
bitrary scale. The F point in the Brillouin zone for phonons 
and the K point in the Brillouin zone for electrons are super¬ 
posed at the origin. Characteristic energies, such as the zone 
width lOeV, maximum phonon energy hwmax 0.2eV, 
Fermi energy ef and energy of acoustic phonons with wave 
vector to kp are specified. Also presented are numerical es¬ 
timates of these quantities. The gray band confines the tem¬ 
perature range of applicability of the theory developed. The 
value O.leV for Fermi energy was chosen because it provides 
the degeneracy of electron gas a room temperature (0.026eV). 
The energy of acoustic phonons with comparable wave vectors 
is 0.002eV, which corresponds to 23 K. It separates equipar- 
tition and Bloch-Griineisen temperature regimes. 


eters are measured in an experiment. Through them can 
be expressed all the results obtained. Doping is another 
way to control the Fermi level in graphene^. 

A temperature gradient in the sample produces a di¬ 
rected phonon flow (the so-called ” phonon wind”). The 
behavior of the phonon system in the presence of a 
temperature gradient can be described in the relaxation 
time approximation with the Boltzmann kinetic equa¬ 
tion. The kinetic equation describing the correction 
= ^phiq) - to the equilibrium distri¬ 

bution function for phonon mode occupation numbers in 
the presence of temperature gradient reads as 


qdKhiQ) 

"""q dT 


kBVT = 
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( 1 ) 


and, hence, the correction itself can be written as 


( 2 ) 

where Tph{q) is the relaxation time of phonon with wave 
vector q. 

Here the relaxation time should be taken for the in 
plane acoustic phonons with a wave vector magnitude 


from 0 to 2kp. Within the temperature range treated 
in this letter, phonon scattering from point defects and 
sample edges are dominant inechanism o^'^d^i^o . For a 
sample 2 fim in size and with the same concentration of 
point defects as observed in Ref»i, scattering from edges 
dominates over that from point defects. The relaxation 
time for edge scattering can be estimated as the ratio of 
the minimum linear size of the sample L (de-facto mean 
free path) to the sound velocity Vs- For the corresponding 
time for a 2^m sample one obtains 100 ps, and it does not 
depend on the magnitude of the phonon wave vector q 
for in plane phonons, thus permitting us to write simply 
in what follows Tph- 

Thermal conductivity of graphene reaches the value of 
5000 Wm“^K“^— In a general case, the relaxation 
time of phonons that generates the drag effect can not 
be straightly connected with thermal conductivity of the 
sample because the phonons from whole Brillouin zone, 
the phonons from other branches and other mechanisms 
of phonon scattering (including the three-phonon Normal 
and Umklapp processes at the temperature higher than 
300 K) define total thermal conductivity of the sample. 

However, the value of the relaxation time of the 
phonons generating the drag effect can be connected 
with low temperature (T < 30 K) thermal conductivity 
of the sample. It is known that the thermal conductiv¬ 
ity of graphene at low temperatures is provided by ZA 
phonons and corresponding relaxation mechanism is scat¬ 
tering on the edges of the sample. The dispersion rela¬ 
tion of ZA phonons is quadratio^: uJzA{q) = Kq^, where 
K = 3.13 • 10“^cm^s“^— . The relation (9) from! after 
substitution the relaxation time TzA{q) = L/vzA{q) and 
some algebra connects the phonon mean free path L and 
thermal conductivity k: 


K = 0.35 


fcV2y3/2 

Kl/2hH3/2 


L, 


( 3 ) 


where h ~ 0.35 nm estimates the thickness of 

graphene layer—. This expression gives approximately 
900 Wm“^K“^ for the sample with L = 2^m at 50K. 

The rate of electron transitions from a state with wave 
vector k into a state with wave vector k' = k + q driven 
by scattering from a phonon with wave vector q (involv¬ 
ing absorption or emission of one acoustic phonon) is 
described by the ” Fermi golden rule”: 
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where M^-ph (?) is the matrix element of electron scatter¬ 
ing from a phonon specifying the magnitude of electron- 
phonon interaction, M is the carbon atom mass, S the 
sample area, and is the number of elementary cells in 
the sample. The combination of the three last quantities 























3 


can be expressed through the surface density of graphene 
p = 6.5- 10“®gcm“^ and the lattice constant ag = 2.46A. 

For long-wavelength acoustic phonons the modulus of 
the matrix element depends on the phonon wave vec¬ 
tor and the scattering angle 0k+q,k as Mf.-ph{(l) = 
Dq^{\ + cos(0fe+q.fc)) ^5,27,37 _ term in brackets sup¬ 
presses electron backscattering^^. In the case of com¬ 
pletely isotropic electron scattering, Me-ph{q) = Dq, the 
phonon drag thermopower would increase fourfold. The 
quantity D can be estimated as 16 eV. Note that this 
value is very close to the combination hvF‘2Trlag, which 
permits the estimation of the graphene electron bands. 


Should a more accurate expression for matrix element 
Me-ph as a function of the wave vectors of electrons and 
phonons be required, it can be readily derived through 
the deformation potential precisely calculated in Refi^. 

In the first order of perturbation theory the contribu¬ 
tion to the drag comes from 4 types of electron tran¬ 
sitions, which correspond to the coming and leaving of 
electrons into the state with wave vector k following ab¬ 
sorption and emission of a phonon with wave vector q. 
The rate of variation of the electron distribution function 
mediated by interaction with phonons can is represented 
in Eq. ([5]) in full form in the standard way through the 
phonon collision integral^^. 


( 
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In Eq. (O (gfe) is the equilibrium distribution func¬ 
tion of electrons with a chemical potential depending on 
ep and temperature. At low temperatures, the chemi¬ 
cal potential can be assumed equal to sp, and the cor¬ 
rection accounting for temperature variation may be ne¬ 
glected. Integration is run over the phonon wave vector 
q over the whole Brillouin zone. Scattering from acous¬ 
tic phonons does not initiate interband and intervalley 
transitions and variation of the electron spin. 

Obviously enough, an equilibrium distribution func¬ 
tion describes randomized phonons, and drag is medi¬ 
ated only by the nonequilibrium term which should 
be substituted from © into 0, and further, into 0. 
vp being approximately 50 times larger than Vs, elec¬ 
trons scatter from phonons in an almost elastic way. The 
energy conservation law is fundamentally important in 
including into (O the distribution function of electrons 
before and after their scattering from a phonon^. One 
can set 

fW^k + q) = f^°Hk)±^^hlOpn{q). ( 6 ) 

The sign identifies here absorption of a quantum of 
sound, and the sign, emission. Conversely, integrat¬ 
ing in angle, one can set 

which corresponds to neglecting phonon energy in a con¬ 
sidering of scattering. As a result, interaction with a 
phonon flux transforms the collision integral to 


where the dimensionless parameter is introduced 


Gph{k) 


S L 

2MNu h^Vp Vs 
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We readily see that the phonon collision integral in¬ 
cludes the quantity and a scalar product of the 

electron wave vector k by the vector describing exter¬ 
nal factor Vr. Hence, the phonon collision integral can 
be interpreted as a field term of the Boltzmann kinetic 
equation. This makes possible calculation of the nonequi¬ 
librium term (fc) for the electron distribution func¬ 
tion in the relaxation time approximation by 

straightforward multiplication of the right-hand part of 
relation relation ([5|) by the electron transport relaxation 
time r(efe). 

One can now calculate the electric current, as well as 
the current j generated by the external electric field. In 
derivation of the current j = ^ J dkvkf^^\k), the first 
approximation (delta-function approximation) in relation 

0 for ^ will be already sufficient. The final expres¬ 
sion for the contribution provided by the phonon drag to 
thermopower looks as 


^ p^kphLI' D y/ ep y 
2eMagVs \hvpTT/ag) \hvpT:/ag) 


The transport electron relaxation time cancels when 
current is divided by the electrical conductivity of the 
sample. 

The most remarkable feature in the relation obtained 
is that it does not contain an explicit dependence of ther¬ 
mopower on temperature, which argues with2^. The ob¬ 
tained relation correlates with the observation that only a 
















4 



0,02 0,04 0,06 0,08 0,10 0,12 0,14 0,16 0,18 0,20 


FIG. 2. Two contributions to the Seebeck coefficient for a 
sample 2 /rm in size and regime of phonon scattering at the 
boundary plotted as a function of Fermi energy. Solid curves 
relate to room temperature, more specifically, the solid black 
line visualizes the phonon drag contribution Sph, and the solid 
grey line, the diffusion contribution Sd- These graphs start 
at O.leV, where the degenerate electron gas approximation 
holds. The dotted black line depicts the phonon drag contri¬ 
bution, and the dotted grey line, the diffusion contribution at 
T = 77K. 

directed phonon flux (phonon wind) contributes to ther¬ 
mopower; now, in accordance with relation 0, this flux 
is proportional to temperature gradient and does not de¬ 
pend directly on the temperature itself. 

The non-dimensional factor in. relation (Uni is 
close in magnitude to 8, and and the ratio ap¬ 

proaches 2. 

The obtained phonon drag thermopower should be 
compared with well known formula (see for instance 
relation (4) from Refi^) of diffusion thermopower in 
graphene: 


In the case of a degenerate electron gas, the diffusion 

thermoelectric current turns out to be zero in the first 
0 ^( 0 ) 

approximation for —gj—. Taking the second term in the 
Sommerfeld expansion results in the proportionality to 
T/ep, which figures in (fTTjl . This expression for the dif¬ 
fusion factor in the thermopower in graphene correlates 
with the Mott formula and coincides with that devel¬ 
oped for metals^ and graphite^. To bring the result in 
compliance with the mechanism of carrier scattering in a 
particular sample, one has naturally to introduce a factor 
of order unity. 

Figure [5] plots the dependencies of the contribution to 
thermopower due to electron drag by phonons and of the 
contribution to thermopower due to diffusion, as func¬ 
tions of the Fermi level, which were calculated by rela¬ 
tions m and (ED, respectively. 


In most of the available relevant publicationsiSiiirJZ 
thermoelectric properties of graphene find straightfor¬ 
ward interpretation in terms of the Mott formula. No 
phonon drag was observed in these studies. Small devia¬ 
tions from the Mott formula were recorded, however, at 
low temperatures. These deviations should most proba¬ 
bly be assigned only to a decrease of measurement accu¬ 
racy. Only the data reported in RefiiS may be interpreted 
as reflecting saturation of thermopower with decreasing 
temperature. Finally, several graphs in RefJ^ suggest 
growth of thermopower with increasing Fermi level. 

The first reason for the phonon drag effect becoming 
suppressed could be the reverse effect exerted by elec¬ 
trons on the phonon system, accompanied by the corre¬ 
sponding decrease of the phonon relaxation time^ The 
second possible factor could be the effect of the sub¬ 
strate, which may likewise reduce the phonon mean free 
paths. Intrinsic phonons of graphene may be expected 
to become scattered from phonons of the substrate or 
from substrate irregularities. For instance, the limita¬ 
tion of the electron mobility caused by the phonons from 
the substrate was indicated in^. Molecular dynamics 
modeling^ yields a value of about 5ps for the relaxation 
time of long-wavelength acoustic phonons in 4 x 4nm^ 
graphene sheets on a Si02 substrate. For a free-standing 
sample, the relaxation time of such phonons turned out 
to be about 10 times larger, approaching in magnitude 
the value of 100 ps used in this study. Nevertheless, it 
appears pertinent to note that interaction of mechani¬ 
cally exfoliated graphene with substrate is assigned to 
weak van-der-Waals forces^i^. Substrate, in particular, 
does not affect the shape of Raman spectra^^, which are 
very sensitive probes of the vibrational characteristics of 
a crystal. 

A calculation was made of the phonon drag contribu¬ 
tion to thermopower in graphene in the Bloch-Griineisen 
(BG) temperature regimei^. The BG regime specihes 
the case where temperature is not high enough to drive 
excitation of phonons with a wave vector on the order 
of the characteristic electron wave vector. The Seebeck 
coefficient in this case is proportional to T^, the situa¬ 
tion coinciding with the case of graphite^^. To realize 
the BG regime in graphene, however, the temperatures 
should be at their lowest limit. As already mentioned, in 
the case of a degenerate electron gas phonons interact¬ 
ing with electrons of energy on the order oi ep ^ O.leV 
have a a temperature of about 20 K. The very noticeable 
predictions of thermopower enhancement are made for 
graphene with band gap opened for both diffusion^i and 
phonon drag contributions^. 

As it was mentioned above, the phonons that generate 
the drag of electrons are very small part of the excited 
phonons that define thermal conductivity. In theory it 
allows to manipulate the thermal conductivity and the 
phonon drag thermopower separately. We suggest that 
the phonon drag thermopower is defined by the scatter¬ 
ing of phonons on the edges of the sample. Thermal 
conductivity can also rise nearly up to the room tern- 
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perature with increasing the sample size. But due to 
relation ZT = aTS^ /n for thermoelectric hgure of merit 
it will grow with increasing of sample size linearly or even 
faster. The problem of the thermal conductivity reduc¬ 
ing is complicated because making samples with point 
defects, placing samples onto the substrate or increasing 
ambient temperature will reduce both electrical and ther¬ 
mal conductivity of the sample and these effects should 
be considered simultaneously. Thereby tuning the size of 
the sample seems to be the straightest way of improving 
graphene thermoelectrical properties. 

The quality of graphene samples prepared is improv¬ 
ing continuously^ii^iiS. As demonstrated in this study, 
in a high quality sample with a large phonon mean free 
path, which probably correlates with high thermal con¬ 
ductivity, at temperatures of the environment from liquid 


nitrogen to room temperature, and at a high Fermi level 
one may expect the contribution to thermopower due to 
phonon drag to dominate over the diffusion-associated 
one, and a sharp increase of the Seebeck coefficient and 
thermoelectric figure of merit. 
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